NT@UW-02-029 



Chiral 1/M 2 corrections to B^*> — > D^*> at Zero Recoil in 
Quenched Chiral Perturbation Theory 

Daniel Arndtf] 

Department of Physics, Box 1560, University of Washington, Seattle, WA 98195-1560, USA 

(Dated: February 7, 2008) 

Abstract 

Heavy quark effective theory can be used to calculate the values of the semileptonic B^ — ► £)W 
decays in the limit that the heavy quark masses are infinite. We calculate the lowest order chiral 
corrections, which are of 0(1/M 2 ), from the breaking of heavy quark symmetry at the zero recoil 
point in quenched chiral perturbation theory. These results will aid in the extrapolation of quenched 
lattice calculations from the light quark masses used on the lattice down to the physical ones. 
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I. INTRODUCTION 



The Cabbibo-Kobayashi-Maskawa (CKM) matrix describes the flavor mixing among the 
quarks, its elements are fundamental input parameters for the standard model. Their precise 
knowledge is not only crucial to determine the standard model but also to shed light on 
the origin of CP violation. The matrix element that parameterizes the amount of mixing 
between the b and c quarks, V^, can be extracted from the exclusive semileptonic B meson 
decays B — > Dlv and B — > D*lv, where I = e, fi. Heavy quark effective theory (HQET) 
(for a recent review, see [|TJ), which is exact in the limit of infinite masses M for the heavy 
quarks, predicts the width of the process B — > D*lv as 

rIV it/ I 2 

£(* - *) = ^KM^ D -M 2 , (i) 

where u — v ■ v' is the scalar product of the 4-velocities v and v' of the D* and B mesons, 
respectively. K(oS) is a known kinematical factor and T(u) is a form factor whose value 
at the kinematical point to — 1 is ^"(1) = 1 in the M — > oo limit. There are, however, 
perturbative and non-perturbative corrections to ^(l), 

Fb^d*(1) = Va + 5 1/M 2 + ..., (2) 

where the parameter rj^ « 0.96 is a QCD radiative correction known to two-loop order 
and Si/m 2 are non-perturbative corrections of 0(1/ M 2 ) to the infinite mass limit of HQET. 
Note that, according to Luke's theorem || there are no 0(1/ M) corrections at zero-recoil. 
One chooses the zero-recoil point because, for u = 1, Tb^d* can be expressed in terms of a 
single form-factor Jia 1 given by 

This is in contrast to the general case u > 1 for which Tb^d* (^) is a linear combination of 
several different form factors of B — > D*lv mediated by vector and axial vector currents. 

Several experiments, most recently by CLEO ||, have determined the product 
(Tb^d*(X)\ Vcb\) 2 by measuring cITb^d* /duj and extrapolating it to the zero-recoil point. 
The mixing parameter \V c b\ can then be extracted once the value ^^^.(l), that encodes 
the strong interaction physics, has been evaluated. The uncertainty in \V c b\ is therefore 
determined by the experimental errors and by theoretical uncertainties in the determination 
of J-'b^d*(1)- Presently, the theoretical uncertainties dominate. 1 

A model-independent way of calculating J-"(l) is provided by numerical lattice QCD sim- 
ulations. In this method one implements field theory non-perturbatively using the Feynman 
path integral approach. The fermion determinant that arises from the path integral is very 
costly to calculate; it is often set to one in an approximation called quenched QCD (QQCD). 
This corresponds to dropping the contribution from virtual quark loops which are made of 



1 Similarly, one can use the decay B — > Dlv to extract (.FB^D(l)|V r c b|) 2 from the measured dTB^o/doj. 
However, cLTb^d / dui is more heavily suppressed by phase-space near to = 1 than dTs^D' I du). In addition, 
the B — > D channel is experimentally more challenging. Thus the extraction of \V c b\ from this channel is 
less precise but serves as a consistency check. 
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"sea" quarks that have propagators not connected to the inserted external operators. "Va- 
lence" quarks, those that are connected to the inserted operators, however, are kept. Re- 
cently, such calculations have been performed || |6|, |7|, |J for the decays B — > D^Hv. Several 
systematic uncertainties, such as from statistics and lattice space dependence, contribute 
to the error of these calculations. Another contribution to the uncertainties comes from 
the chiral extrapolation of the light quark mass. Since lattice QCD simulations are limited 
by the available computing power they presently cannot be performed with the physical 
masses of the light quarks. Therefore one needs to extrapolate from the heavier masses used 
on the lattice (of order the strange quark mass) down to the physical light quark masses. 
This extrapolation can be done by matching QQCD to quenched chiral perturbation theory 
(QxPT) and calculating the non-analytic corrections 5um 2 in Eq. (0) in Q%PT. The for- 
mally dominant contributions to these corrections come from the hyperfme mass splitting 
between the heavy pseudoscalar and vector mesons that stems from the inclusion of heavy 
quark symmetry breaking operators of 0(1/ M) in the Lagrangian. 

In QCD, the corrections due to D meson hyperfme splitting have been calculated in 
chiral perturbation theory (%PT) by Randall and Wise |J. A more complete treatment, 
involving additional corrections due to B meson hyperfme splitting, 0(1/ M) axial coupling 
corrections, and 0(1/ M) corrections to the current, has been given in [|Kj]. Recently, the D 
meson hyperfme splitting corrections have also been determined in partially quenched chiral 
perturbation theory (PQ^PT) [pT[] for partially quenched QCD (PQQCD). In contrast to 
QQCD, PQQCD does not drop the contributions from sea quark loops but gives different 
(and separately varied) masses to the sea and valence quarks. Usually, sea quarks are heavier 
than valence quarks, so that the fermion determinant — although no longer equal to one - 
is much less costly to calculate than in ordinary QCD. 

In this paper we calculate the 0(1/ M 2 ) corrections in Eq. (|2]) due to D and B me- 
son hyperfme splitting in quenched chiral perturbation theory (Q;\PT). As we will show, 
these corrections are — upon expanding in powers of the hyperfme splitting A — of order 
Aqcd 3 ™^ 2 / (M n m^ 2 ) for n > 2 and formally larger than those coming from the inclusion of 
0(1/M) heavy quark symmetry breaking operators in the Lagrangian and current which 
are suppressed by powers of Aqc D /M. This argument is similar to the one that applies to 
XPT0. 

Our QxPT calculation can be used to extrapolate lattice results M that use the quenched 
approximation down to the physical light quark masses. So far, this extrapolation has been 
based upon the ;\RT calculation |5) . Using QxPT should therefore give a better estimate of 
the uncertainties related to the chiral extrapolation. 

A central role in the lattice calculation of B — > D* |7|, § is played by the double-ratios 
of matrix elements 

(D\c 1 %\B)(B\h°c\D) 
+ (D\c^c\D)(B\h%\BY 1 ' 

= (D*\cy°b\B*)(B*\h°c\D*) 
1 {D*\afc\D*)(B*\b7%\B*Y [ ' 

and 

n = (D*\nh 5 b\B)(B*\h j l5c\D) () 
Al {D*\&y^ 5 c\D)(B*\bj^ 5 b\B)' 1 ' 

Since the numerator and denominator are so similar, statistical fluctuations are highly corre- 
lated and cancel in the ratios to a large degree. The 0(1/ M 2 ) correction to the double ratios 
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can therefore be calculated fairly accurately and used to derive the 0(1/ M 2 ) correction to 
the matrix elements themselves. For this reason, we also calculate 0(1/M 2 ) corrections 
to the decay B* — > D* in addition to the experimentally accessible decays B — > D and 
B — > D*, and thus the corrections to 1Z + , IZi, and TZa 1 - 



II. QUENCHED CHIRAL PERTURBATION THEORY 

In a world where all the quark masses are large compared to Aqcd internal quark loops 
are suppressed and the results from QQCD are close to those from QCD. In the real world, 
however, light quarks are light (<^C Aq CD ) and contributions from internal quark loops are 
substantial. One can nevertheless study the low-energy behavior of QQCD by its effective 
low energy theory, Q^PT [0, |E| [TJ, [LJ, |T§ . 



We consider a theory constructed from three light valence-quarks, u, d, s, and three light 
bosonic quarks, u, d, s governed by the Lagrangian of QQCD 

C= q a W-rn g y a q a + £ f(i$ - m^)\q~ a = £ Q\iJp - m Q )]Q j . (7) 

a=u,d,s a=u,d,s j=u,d,s,u,d,s 

Here both types of quarks have been accommodated in the six-component vector Q with 
the three quarks q a in the upper three entries and the three bosonic ghost-quarks q~a in the 
lower three entries. The graded equal-time commutation relation governing the valence- and 
ghost-quarks is 

QU*)Qf (y) - (-i)^'Qf (y)0f (x) = rt^ 3 (x - y ), (8) 

where a and f3 are spin- and % and j are flavor-indices. The graded equal-time commutation 
relations for two Q's and two Q^s are analogous, rjk is given by 

fl fork =1,2,3 
lk \ for fc = 4, 5, 6 ' 1 ' 

The quark mass matrix is given by mq = diag(m M , m^, m s , m u , m^, m s ), i.e., the fermionic 
and bosonic quarks have equal masses but different statistics. Therefore the contributions 
of fermionic and bosonic quarks in virtual quark loops cancel exactly. The Lagrangian in 
Eq. (0) exhibits a graded symmetry [S77(3|3)l ® S77(3|3).r] x U(l)y that is assumed to 
be broken spontaneously to SU(3\3)v x U(l)v- The dynamics of the emerging 36 pseudo- 
Goldstone mesons can be described at lowest order in the chiral expansion by the Lagrangian 

£ = L-sti (d^d^E) + A str (m Q £ + m^S f ) + ad^ d^ - ^ (10) 

where 



and 

X 7T 



$ = ( 71 X ~ ) ■ (12) 
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Here the tt, tt, and \ are 3x3 matrices of pseudo-Goldstone bosons with quantum numbers of 
qq pairs, pseudo-Goldstone bosons with quantum numbers of qq pairs, and pseudo-Goldstone 
fermions with quantum numbers of qq pairs, respectively, 



TT 



T] u 7T + K + 
71- T]_ d K° 

K~ K° r] s 



TT 



\k- k° fj s J 



and x 



Xvu Xn+ Xk+ 

Xvr- Xfld Xk° 

Xk- Xk° Xvs 



(13) 



The pion decay constant is fixed by Eq. (|Tl"D and / = 132 MeV in QCD. 
The flavor-singlet field $ is defined as 

1 l 

$0 = 7E stT ^ = 7f (r/, ~^ 



(14) 



where str() denotes a supertrace over the flavor indices. <3>o is invariant under [S77(3|3)z, <g) 
SU(3\3)r] x U(l)v and thus arbitrary functions of it can be included in the Lagrangian. 
To lowest order in the chiral expansion only the two operators included in Eq. fllOf) with 
parameters a and fi remain and are understood to be inserted perturbatively flip . Notice 
that this singlet field $o is n °t heavy as in %PT and therefore cannot be integrated out. It 
introduces a new vertex, the so-called hairpin. 

Upon expanding the Lagrangian in Eq. (|T0|) one finds that the mesons with quark content 
mm, dd, and ss, the only ones relevant for our calculation, have masses given by 



m 



qq 



8\m q 



(15) 



III. INCLUSION OF HEAVY QUARKS 

The .D-mesons with quantum numbers of cQ can be written as a six-component vector 

D = (D w D d ,D„D a ,D i ,D s ). (16) 

Heavy quark symmetry is provided by combining creation and annihilation operators for the 
pseudoscalar and vector mesons, D and D* respectively, together into the field H D 



D 



H 



D 



7 



H 



1 + 



(17) 
(18) 



where v denotes the velocity of a heavy meson. In HQET the momentum of a heavy quark is 
only changed by a small residual momentum of O(Aqcd)- Hence, v is not changed and H is 
usually denoted by an index v which we have dropped here to unclutter the formalism. In the 



heavy quark limit, the dynamics of the heavy mesons are described by the Lagrangian [17, 18 



C 



D 



tv[H^v^5 ab + iV£)H?\ + g tr(H°H° lul5 A» ba ) + 7 tr(H?H?w s ) BbrA" (19) 
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where the traces tr() are over Dirac indices and supertraces str() over the flavor indices are 
implicit. The additional coupling term involving $o ~ strA M is a feature of Qx?T and not 
present in xPT. The light-meson fields are 

4. = - e^e f ) = ~- f o^ + 3 ) (20) 

and 

^ = ^(^ + #U f ) = 2^[$, 5^] + 0($ 4 ). (21) 

Expanding the Lagrangian Co to lowest order in the meson fields leads to the (derivative) 
couplings DD*d(p and D*D*d<p whose coupling constants are equal as a consequence of 
heavy quark spin symmetry. The DDd<p coupling vanishes by parity. 

An analogous formalism applies to the fields B and B* which are combined into H B . 
Note that the axial coupling g is the same for H D and H B mesons as dictated by heavy 
quark flavor symmetry. 

We do not include terms of order m q ~ y/m n in the Lagrangian as explicit chiral symmetry 
breaking effects are suppressed compared to the leading corrections. The presence of these 
terms is implied by the nonzero masses m qq . 

IV. MATRIX ELEMENTS OF flW _» D^W 

The non-zero hadronic matrix elements for B^*' — > can be defined in terms of the 16 
independent form factors h±, hy, hA 123 , and /11...10 as |1|, |T^ 



(mmm = h+M{v + vY + Mw)( „ _ vYi (22) 

yJm B m D 

(j^mm = -m«)— «x«„ (23) 

y/m B m D * 

(D%v',e ')\cri 5 b \B(jj))_ = _ iflAi ^ UJ + 1 y^ + ih ^ v . e '*) v ^ + ih M (u){ve'*)v'^ (24) 
^/m B m D * 



{D\v\e ')\cYb\B 2i^ri = _(^. e )[ Al ( w )( v + v y + A2 ( w ) (v _ v y] + ^ 3 ( w )( 6 «.. v )^ 
^/m B *m D * 

+h 4 {uj){t ■ v')e'^ - (e • </)(e'* • u)[/i 5 (u;K + /^(a;)^], (25) 



^ aK<5 {tlt 5 [h 7 (u){v + v'Y + h 8 (u){v - i/)"] 



and 

(D*(v', t')\c^^b\B* (v , t)) 
y/m B *m D * 

+ t/V [h 9 (uj)(t'* ■ vy + h w (u)(t ■ v')t'^}} .(26) 

Here, u — v ■ v' and v [t] and v' [t'\ are the velocities [polarization vectors] of the initial state 
B^ meson and final state meson, respectively. Note that we will not explicitly calculate 
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matrix elements of B* — > D as these can be easily related to the B — > D* calculation by a 
hermitian conjugation of the matrix elements and an interchange of the c and b quarks, i.e., 

£(*) «_> £)(*). 

In the heavy quark limit the matrix elements in Eqs. (p2|) (p6|) are reproduced by the 
operator 

af{l - l5 )b -> -C(o;)tr[^V(l - 7*)^*]- (27) 

Here, £(co>) is the universal Isgur-Wise function |p0[ , [21|| with the normalization £(1) = 1. To 
lowest order in the heavy quark expansion one finds 

h + (uj) = h v (uj) = h Al (uj) = h Ai (uj) = hx(uj) = h 3 (co) = h A (uj) = h 7 (uj) = £(lu) (28) 

and the remaining 8 form factors vanish. 

The discussion of the B^ — > D^W matrix elements is similar for different flavors of the 
light quark q content of the B^ and mesons; it applies equally to q = u, d, or s as 
the theory splits into three similar copies of a one-flavor theory. In the limit of light quark 
SU(3)v flavor symmetry the matrix elements (and in particular the Isgur-Wise function) 
are therefore independent of the light quark flavor. However, in nature the masses of the u, 
d, and s quarks are different and 577(3) y is not an exact symmetry. Therefore our results 
will include terms that depend upon m q via the meson masses m qq defined in Eq. (ITS). 



V. l/M CORRECTIONS 



The lowest order heavy quark symmetry violating operator that can be included in the 

H n D a^H?a h 



Lagrangian in Eq. (fi9| ) is the dimension-three operator ^-tr 



It vio- 



lates heavy-quark spin and flavor symmetries and comes from the QCD magnetic moment 
operator ca fJ,u G^ lJ T A c, where G A is the gluon field strength tensor and T A with A = 1 ... 8 
are the eight color SU(3) generators. This operator gives rise to a mass difference between 
the D and D* mesons of 

A D = m D , - mD = -8-J^. (29) 
This effect can be taken into account by modifying the D and D* propagators which become 



and ™ao\wv f30) 
2(v k + 3A D /4 + ie) 2(v ■ k - A D /4 + it) ' V ; 

respectively, so that in the rest frame, where v = (1,0,0,0), an on-shell D has residual 
energy of SA^/A and an on-shell D* has residual energy of Ad/ 4. A similar effect due 
to the inclusion of a QCD magnetic moment operator for the b quark applies to the B^ 
mesons. 

There are no corrections to the matrix elements for the semileptonic decays B^ — > D^*'ev 
of 0(1/ M) at zero-recoil according to Luke's theorem M. The leading corrections enter at 
0(1/M 2 ). In addition to tree-level contributions from the insertion of 0(1/M 2 ) suppressed 
operators into the heavy quark Lagrangian or the current there are one-loop contributions 
from wavefunction renormalization and vertex correction. These one-loop diagrams have a 
non-analytic dependence on the meson mass m qq and depend on the subtraction point \i. 
This dependence on \i is canceled by the tree-level contribution of the 0{1/M 2 ) operators. 
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(b) 

FIG. 1: Graphs contributing to wavefunction renormalization for heavy (a) pseudoscalar and 
(b) vector mesons. A thin [thick] line denotes a heavy pseudoscalar [vector] meson, a dashed line 
denotes the <&o, while a dashed-crossed line denotes the insertion of a hairpin. A full [empty] vertex 
denotes a g [7] coupling. 



Because of the absence of disconnected quark loops in QQCD, which manifests itself as a 
cancellation between intermediate pseudo-Goldstone bosons and pseudo-Goldstone fermions 
in loops in QxPT, the only loop diagrams that survive are those that contain a hairpin 
interaction or a 7 coupling. 

The wavefunction renormalization contributions for the pseudoscalar and vector meson, 
Zjj/b and Zp/ B , respectively, come from the one-loop diagrams shown in Fig. ([]]) and have 
been calculated in |T7L pl| p2fl. Including the a coupling we find for these diagrams 



l + i^AJ-^AJ+^CA) (31) 



and 



!gi ffl( _ A) _^ 2( _ A)+ ^, 



+^.(0)-^,(0) + ^F I (0). (32) 

The functions H\, H2, and F\ come from loop integrals and are given in the Appendix. Note 
that in the heavy quark limit where A = one recovers Z = Z*, as required by heavy quark 
symmetry. 

The vertex corrections come from one-loop diagrams. The non-vanishing contributions 
are shown in Fig. @. Combining the wavefunction renormalization and vertex corrections 
and including a local counterterm to cancel the dependence on the renormalization scale \x 
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B B* D* D 



i — ^— 

B B* D* D 



■x- N 



(a) 



B B* D D* 



B B* D D* B B* D* D* B B* D* D* 



(b) 



B* B D D* B* B D D* B* B* D* D* B* B* D* D* 



FIG. 2: QxPT graphs which contribute to the vertex correction of the form factors (a) h + (l), (b) 
/^(l), and (c) hi(l). A full [empty] square denotes the insertion of the operator 07^756 [cj^b]. 



we find the following corrections for the form factors: 
5h + (l) 



-MaO + o + n 



ff 5 (A B , A D ) - a tf 8 (A B , A D )} - ^G 5 (A B , A D ) 



6^7 



xm + 



1 



(4vr/) 2 V3m 2 
+0({A i? ,A D } 3 ), 



# 2 a 



171 , , 

Iog^)(A 



b-A d ) 2 



(33) 



^(1) 



,2 



1 



+ 



2<r 



3/ 2 
2ig* 



2 2 
2 if 5 (A B ,-A D )-aif 8 (A B ,-A D 



^G 5 (A B ,-A D ) 



H 5 (A B ,0) -aH 8 (A B ,0)} - -jLLG 5 (A B ,0) 



3/ 2 

XM+ (A7cfy \9m* 
+0({A B ,A D } 3 ), 



v g 2 g _ 2gi_ 
9 3 



log — (3A| + A 2 D + 2A B A D ) 



(34) 
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and 



5h(l) = x^) + £*_J: + £*1_L 



/(r r.-a^A^-A^-a^C-AB.-Afl)] -^G 5 (-A B ,-A D ) 



2i 5 2 r 2 „ / n tr /n n\l 4i #7 



tf 5 (0, 0) - a tf 8 (0, 0)] - -^G 5 (0, 0) 



^ ^2 |_' U O \ 5 / o \ 7 /J J? 

Xi(/i)+ - V« 2, 7 



1 (qM 

(4tt/) 2 V9m 2 
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log^)(A B -A D ) 2 
A 4 



+G({A B ,A D } 3 ), (35) 

which are defined by h + (l) = l + 5h + (l) and analog expressions for 5/^(1) and The 
functions if 5, H 8 , and G5 come from loop-integrals that are listed in the Appendix and we 
have defined m = m qq . The insertions of tree-level 0(1/M 2 ) operators are represented by 
the functions Xa x (h), and Xi(fj) which are independent of m and exactly cancel the 

jj, dependence of the logarithm. These functions can be extracted from lattice simulations 
by measuring the zero-recoil form factors for a varying mass of the light quark. 

Experimentally, Ad pa 142 MeV and A B pa 46MeV so that the ratios Ap/m and A B /m, 
which enter the form factor corrections through the function R(A/m) (defined in the Ap- 
pendix), are 0(1). On the lattice, however, one can vary all quark masses. Expanding first 
in powers of A and then taking the chiral limit m — > one finds the formal limits given in 



Eqs. (33)-fl35"D where we have only kept the pieces non-analytic in m. This demonstrates 
that the terms linear in Ad and A B , although present in wavefunction renormalization and 
vertex corrections, cancel as required by Luke's theorem ||. The leading order corrections 
areO({A B ,AD} 2 ). 

As a consistency check one can restore heavy quark flavor symmetry by taking A B = Ad- 
Since the 0{1/M 2 ) corrections to h + (l) and h\(l) are proportional to (A B — A D ) 2 they 
disappear as they should since the charge associated with the operators cy^c and bj^b is 
conserved. This argument does not apply for the B — > D* transition matrix element in the 
limit Ad = Ad since there is no conserved axial charge associated with the operators 07^750 
and 67^75 b. 

In the chiral limit, the term proportional to /i 2 , has a 1/m 2 singularity and dominates over 
the terms proportional to a and 7 that are only log- divergent. This is analogous to a term of 
the form {m 2 qq — m 2 j)/m qq found by Savage |lTJ for PQxPT (here, m qq and TTljj £1X6 valence 



and sea quark masses, respectively). In the limit rrijj — > m qq this term, however, vanishes 
as PQxPT goes to %PT where the dominant term is \ogm qq . In Q%PT, on the other hand, 
the l/m 2 qq pole persists, revealing the sickness of QQCD where the hairpin interactions give 
a completely different chiral behavior than in QCD. 

The size of /xq can be estimated from the rj-rj' mass splitting |i~3 |, large Nc arguments |23|, 



4| (Nc being the number of colors), or lattice calculations. These estimates imply /jq 
500 — 900 MeV; for the purpose of dimensional analysis we use //q ~ C(Aqcd)- Taking g ~ 
0(1) we therefore find that 6h+, Sh Al , and 5h x are of order A n /m n ~ A QCD 3n/2 /(M"m™ /2 ) 
for n > 2 and thus larger than tree-level heavy quark symmetry breaking operators that are 
suppressed by Aqcb/M. 

To show the dependence of the zero-recoil form-factors on the mass of the light spectator 
quark it is necessary to know the numerical values of the parameters /Jo, g, a, and 7. In 
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m /m 

q s 



FIG. 3: Dependence of on the mass m q of the light spectator quark in QxPT. For compar- 

ison, the xPT result from B is also shown (dashed line). The result has been normalized to unity 
for m q = m s . We have chosen fiQ = 700 MeV and g 2 = 0.4. 



determining reasonable values for these couplings we follow the discussion by Sharpe and 
Zhang ||18|| . Assuming that g is similar to the %PT value we use g 2 = 0.4. The hairpin 
coupling a is proportional to 1/Nc, and thus assumed to be small; we use two values, a = 
and a = 0.7. The coupling 7 is known to be suppressed by 1/Nc compared to g, the sign is 



undetermined. We take —g < 7 < g (see [|18j and references therein) 



With these parameters, the dependence of h + (l) and /iAi(l) on the mass of the light 
spectator quark m q is show in Figs. (|) and ([|), respectively. The graphs are plotted against 
m q in units of the strange quark mass m s with m q /m s = m 2 /m 2 s where m 2 s = 2m 2 K . The 
behavior of h + (l) in QxPT is dominated at small m by the 1/m 2 pole that is non-existent 
in xPT. Lattice calculations of h + (l) || show a small downward trend for decreasing m q 
down to the chiral limit that is similar to the downward trend seen from the calculation 
(dashed line). The same behavior (down to m q w 0.1m s ) can also be seen for Q%PT for 
a certain choice of parameters (e.g., 7 positiv). The case of /iAi(l) is different as there is 
a pole at m = Afl which is close to the physical pion mass. Here, both D* and 7r can be 
on-shell and the decay B — > D*n becomes kinematically allowed. Lattice calculations of 
^Ai(l) [0 for m q = (0.6 . . . l)m s show a small downward trend for decreasing m q similar to 
the downward trend seen from the %PT calculation [dashed line in Fig. (f3J)]. A similar trend 
down to m q « 0.2m s can also be seen in the Q%PT calculation for a relatively large positiv 
value of 7. 

Although the downward trend in the lattice data for the two cases seems significant as 
the statistical errors are highly correlated, the uncertainty is still relatively high (typically 
±0.01) and the existing lattice data can be accommodated by a wide range of values for the 
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QXPT, a=0.7 
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m /m 

q s 



FIG. 4: Dependence of h Al (1) on the mass of the light spectator quark in Q^PT. The dashed line 
denotes the %PT result H. The numerical values for the parameters are those used in Fig. (p|). 



parameters in the Qx?T Lagrangian. 

Finally, we calculate the double ratios denned in Eqs. 
(El)-©. We find 

K + = l + 2Sh+{l), 



using the results in Eqs. 

(36) 



Tlx = l + 25hi(l) 



(37) 



and 



3f 



{f4 [H 5 (A B , -A D ) + H 5 (A D , -A B ) - H 5 (A D , -A D ) - H 5 (A B , -A B )] 
- a [H 8 (A B , -A D ) + H 8 (A D , -A B ) - H S (A D , -A D ) - H 8 (A B , -A B )}} 



2^7 



, [G 5 (A B , -An) + G 5 (A D , -A B ) - G 5 (A D , -A D ) - G 5 (A B} -A B )\ 



f 



(4tt/) 2 V 9m 2 
+0({A B ,A D } 3 ), 



2g 2 & _ 4^7 
9 3 



log^)(A B -A D ) 2 



(38) 



where X Al (/x) is the counter term associated with lZ Al ■ 
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VI. CONCLUSIONS 



Knowledge of the fiW — > form factors at the zero-recoil point is crucial to extract the 
value of V c b from experiment. In the limit that the heavy quarks are infinitely heavy HQET 
predicts that the form factors h + , h^, and hi are equal, h+(l) = ^Ai(l) = hi(l) = £(!)• 
The formally dominant correction due to breaking of heavy quark symmetry comes from the 
inclusion of a 0(1/M) dimension-three operator in the Lagrangian that leads to hyperfine- 
splitting between the heavy pseudoscalar and vector mesons. These leading order corrections 
are 0({Ab, A^} 2 as required by Luke's theorem. 

Recent lattice simulations using the quenched approximation of QCD have made a big 
step forward in determining these zero-recoil form factors. Presently, however, the simu- 
lations use light quark masses that are much heavier than the physical ones and therefore 
rely on a chiral extrapolation down to the physical quark masses. In this paper we have 
calculated the dominant corrections to the form factors h + , h^, and hi in Q%PT and de- 
termined the non-analytic dependence on the light quark masses via the light meson masses 
m qq . Using these results, instead of the yPT calculation, to extrapolate the QQCD lattice 
measurements of these form factors down to the physical pion mass should give a more 
reliable estimate of the errors associated with the chiral extrapolation. 

We have also calculated the corrections to certain double ratios that are used in lattice 
QCD calculations of the decay B — ► D* . 
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APPENDIX A: INTEGRALS 



We list the functions Hi, H 2 , Fi, H 5 , H 8 , and G$ even though some of them have appeared 
in the literature before UlO , 25fl . Here, m = m qq is the mass of the qq light meson in the loop 
where q = u, d, or s is the light (spectator) quark content of the heavy mesons. We have 
used dimensional regularization with the MS scheme, where 1/e' = 1/e — 7^ + log47r + 1. 
At the end we set 1/e' = 0. As a shorthand we have defined the function 



R(x) 



Hog 



x 



- Vx 2 - 1 



16 



+ \/ x 1 — 1 + ie 



x 



which occurs frequently. We also need its derivative dR/dx given by 

x 



R'(x) 



-R(x)-2. 



(Al) 



(A2) 



x 2 — 1 

For the calculation of the wavefunction renormalization contribution we need the deriva- 
tives of the loop integrals for the diagrams in Fig. ([[]): 

t 



Hi(A) 



m" 1 
log — - - - 1 



4- 



(A3) 
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H 2 (A) 



16tt 2 



A 2 - —m 2 + 2 m 2 - A 2 log — 
3 3 \ fi 2 



+ - Ami? 



m 



(A4) 



and 



F X (A) = 



16tt 2 



10 



f A2 - i™ 2 + < m2 - 2A2 > i log 7? - 7 j + i Amfl I £ 



A 



+ |(A 2 -m 2 )i?' 
3 \m 



(A5) 



For the loop-integrals of the vertex corrections one finds 

^ 1 



H 5 (A,A) = 



m 



16tt 2 



m 

]0g l^-7- 1 A-A 



«| A 

m 



ft i A 

m 



(A6) 



#s(A,A) 



16tt 2 



2m 2 --(A 2 + AA + A 2 ) 
3 



m 



16 



10 



m(5m 2 - 2A 2 ) / A \ m(5m 2 - 2A 2 



m' H = -R 

3 3(A - A) 



3(A - A) 



-R 



A 



m 



,(A7) 



and 



G 5 (A,A) = 



16 ,n> 2 + AA + A2) -5 m2 + 



m 2 --(A 2 + AA + A 2 ) 



m 2 1 

10 g ^~7> 



2m(A 2 -m 2 ) /A\ 2m(A 2 - m 2 ) (A 

H R — R 

3(A-A) \mj 3(A-A) \m 



(A8) 
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